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a'. 

■ We calculate the (parity-violating) spin rotation angle of a polarized neutron beam 
CN . through hydrogen and deuterium targets, using pionless effective field theory up to 
' _h ' next-to-leading order. Our result is part of a program to obtain the five leading 
,_! H ■ independent low-energy parameters that characterize hadronic parity-violation from 

o; 

few-body observables in one systematic and consistent framework. The two spin- 
rotation angles provide independent constraints on these parameters. Our result 
for np spin rotation is ± = [4.5 ± 0.5]radMeV~5 (^g^ 1 '^ + g^~ lp ^ 

[18.5 ± 1.9] radMeV^ (^A/Iof^ ~ 2 9(£zl)^\ while for nd s P in rotation we ob- 
tain ±^jp = [8.0±0.8]radMeV"5 g^- lp ^) + [17.0 ± 1.7] radMeV - * g( 3 ^- 3p ^ + 
[2.3 ± 0.5]radMeV~5 hg { ^°~^ o) - 2g { ^°~^ o) \ where the g {x ~ Y) , in units of 

y—i ■ MeV~2 ; are the presently unknown parameters in the leading-order parity- violating 

j> ■ 

. i-H . Lagrangian. Using naive dimensional analysis to estimate the typical size of the 

i ' couplings, we expect the signal for standard target densities to be -7/p 

[l0~ 7 • • • 10~ 6 ] — for both hydrogen and deuterium targets. We find no indica- 
tion that the nd observable is enhanced compared to the np one. All results are 
properly renormalized. An estimate of the numerical and systematic uncertainties of 
our calculations indicates excellent convergence. An appendix contains the relevant 
partial-wave projectors of the three-nucleon system. 
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I. INTRODUCTION 

Parity-violating nucleon-nucleon interactions cause the spin of transversely polarized 
neutrons to undergo a rotation when traveling through a target medium, even in the absence 
of magnetic fields. In this paper we report the results of a calculation of neutron spin rotation 
from neutron-proton and neutron-deuteron forward scattering using pionless effective field 
theory, EFT(^). This work is part of an effort to provide one consistent EFT(^) framework 
with reliable theoretical uncertainties to parity-violating (PV) interactions in few-nucleon 
systems. We hope that the results presented here, along with the results from Refs. [ll, |2|, 
will assist in the planning, analysis, and interpretation of related PV experiments. 

The PV component of the force between nucleons stems from the weak interactions 
between the standard model constituents of the nucleons. Compared to the parity-conserving 
(PC) part, it is typically suppressed by a factor of ~ 1CT 7 to 10~ 6 ; for reviews see e.g. Refs. [3|, 
4J. Parity- violating neutron spin rotation observables were first discussed by Michel in 
1964 [5j and studied further in Refs. Recently, an upper bound on the effect in 

Helium-4 was obtained at NIST j^. 

At present, the effects of hadronic parity violation on the nuclear level cannot be pre- 
dicted from first principles. Early approaches to PV nucleon-nucleon (NN) interactions 



include the parameterization in terms of S-P wave transitions 



10] and the more widely 



used phenomenological meson-exchange models, particularly the framework developed by 



Desplanques, Donoghue and Holstein (DDH) in Ref. [Hj. The DDH approach provides es- 
timated ranges for PV meson-nucleon couplings based on a number of model assumptions. 
Using the "best values," it has been applied to the study of neutron spin rotation on various 



targets 
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12Hl5|. Calculations have also been performed in the so-called hybrid formalism 
17] . where phenomenological wave functions in the strong sector are combined with a 
PV effective field theory (EFT) treatment. We discuss this further in Sec. IVII1 

In order to avoid model assumptions and to treat all interactions within a unified frame- 
work, we apply effective field theory techniques consistently to the neutron-hydrogen and 
neutron-deuterium systems. The identification of a small parameter allows a systematic ex- 
pansion of our results and a reliable estimate of the size of theoretical errors. In particular, 
since typical neutron energies in parity-violating spin-rotation experiments on the lightest 
nuclei are low enough that pion exchange cannot be resolved, we use the pionless EFT with 



3 



only nucleons as dynamical degrees of freedom. This theory has proven highly successful 



in the parity-conserving sector; for an overview see e.g. Refs. 
scriptions of hadronic parity violation see Refs. 121 
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-|20|. For early EFT de- 



4231 ] . A comprehensive formulation of P V 



effects in EFTs with and without pions was given in Ref. 24]. 

In the PV sector of EFT(^), five independent operators appear at leading order (LO) 
in EFT(^f). They correspond to the five transition amplitudes from Ref. 10J expressed 
in a field theory language. The five accompanying parameters, or low-energy constants 
(LECs), encode the unresolved short-distance physics. At present, only experimental input 
can determine these couplings without introducing additional model dependence. Note, 
however, that for the pionful sector a first study using lattice simulations to determine 
the PV ttNN coupling has been performed 25|]. Our final, next-to-leading order (NLO) 
results for np and nd spin rotation provide these processes in terms of the PV LECs, along 
with estimates of the associated theoretical errors. Measurements of these observables can 
determine two independent combinations of the PV low-energy constants. 

This article is organized as follows: We first review the general formalism of neutron 
spin rotation in Sec. [Ill and present the necessary PC and PV pieces of the Lagrangian in 
Sees. II III and IIV| respectively. The results for neutron spin-rotation on the proton up to 
NLO, along with error estimates, are given in Sec. |V] The results for deuterium up to NLO 
are derived with a detailed discussion of the expected theoretical uncertainties in Sec. I VI I In 
Sec. I VI 1 1 we estimate numerical predictions and compare with earlier work. Conclusions and 
outlook are given in Sec. IVIIII Appendices contain the general construction principle for the 
partial-wave projectors of the three-nucleon system and its results for the S- and P-waves 
as well as details of the numerical calculations. 1 



II. NEUTRON SPIN ROTATION GENERAL FORMALISM 



In this section, we define the spin rotation angle, its relation to the scattering amplitude, 



and the associated conventions we will use. Important resources are Refs. 
that conventions vary. 



27| . but note 



A beam of very-low energy neutrons passing through a medium picks up a phase factor 



1 After completion of this article a paper appeared by Vanasse 26[ , with an analysis of neutron spin rotation 
off the dcuteron using pionlcss EFT to leading order. 
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from scattering in the target. This phase factor is related to the index of refraction n of the 
medium. In the simplest case of very low-energy scattering with plane waves describing the 
incoming state, the phase accumulated after traversing a target of thickness I is given by 

<p = Re(n - l)fci ab Z, (II. 1) 

where is the magnitude of the wave vector of the incoming particle. The index of 
refraction is in turn related to the scattering length a by 

n-l = --^, (11.2) 

where p is the density of scattering centers in the target. In our convention, the scattering 
length and scattering amplitude at zero energy are related by 

M = -—a, (II.3) 

with p the reduced mass of the beam-target system. The phase p can therefore be written 

as 

(p = pl±-R e {M) . (II.4) 

For a beam chosen in the +z direction, a perpendicular polarization in the +x direction 
can be written as 

|x + > = -^(|+> + |-» f (IL5) 

where the states |±) represent states with positive /negative helicity along +z. When tra- 
versing a medium, each helicity state evolves with a phase factor: 

i= (e-^|+) + e-^-R) = (|+> + e^-^R) . (II.6) 

As long as <fi + = 0_, which is the case for parity-conserving interactions, the polarization 
of the beam is unchanged; the state simply picks up an overall phase factor. In the case of 
parity violation, however, <ft + ^ 0_ and the neutron spin is rotated by an amount 

0pv = 0+-0_. (IL7) 

A positive value of the spin rotation angle 0py corresponds to a rotation about the neutron 
momentum in the sense of a right-handed screw. Using Eq. ( 111.4p . the spin rotation angle 
per unit length I is 

1 #PV " Re(M + -M„), (11.8) 



P dl fci a b 

where M.± is the scattering amplitude for ±-helicity neutrons including the statistical mix 



ture of available target spins. For further details see Refs. 
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III. PARITY-CONSERVING LAGRANGIANS AND AMPLITUDES 



A. Two-Nucleon Sector 



A description of EFT(^) and its power counting can be found in Refs. 18M20| , for example. 
Pionless EFT is applicable to energies E < m^/M, where m n and M are the pion and 
nucleon mass, respectively. The short-distance details of the interactions are encoded in the 
low-energy couplings (LECs) of NN contact interactions. Operators and observables are 
organized in terms of the small dimensionless expansion parameter Q = where p typ is 
the typical external momentum or momentum transfer and ~ is the breakdown scale 
of the theory (the scale at which pion exchange can be resolved). The expansion parameter 
is typically 1/5 to 1/3. The following leading-order (LO) and next-to-leading order (NLO) 
calculations only require S-wave interactions in the parity-conserving NN sector; higher 
partial waves are suppressed according to the power counting. It is convenient (in particular 



for three-body calculations) to introduce spin-triplet and spin-singlet dibaryon fie 



d s with the quantum numbers of the corresponding two- nucleon S-wave states 28M30| . The 



ds d t and 



dibaryon (auxiliary) field dt also serves as the deuteron interpolating field since both have 
identical quantum numbers. The relevant terms of the Lagrangian up to NLO are 



C PC =iV f i«9 + 



N-y 




d^(N T P l t N) + H.c. - y [df(N T P s A N) + H.c] 



AM + M 



A, 



CQs 



ldo + m + M 



l 2 s 



+ ... 



using the conventions of Ref. 



{111.1) 



31] . N — ( p ) is the isospin doublet of nucleon Weyl spinors 



p (proton) and n (neutron). With <jj (r^) the Pauli matrices in spin (isospin) space, the 
matrices P t l = ■^-r^a,; and Pf 
1 Sq partial waves [32 ]. 



-j^TiTaVi project two-nucleon states onto the 3 Si and 



The parameters of the Lagrangian are fixed using Z-parameterization 3lJ, [33[ . Choosing 



2 47T 

V =M 



(III.2) 



the LO parameters A s u are determined from the poles of the NN S-wave scattering amp- 
litudes at ry a /t- in the 3 Si-wave, this reproduces the experimental binding energy of the 
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deuteron = — 7 f 2 /M. The leading-order dibaryon propagators are then given by 

Dfifo, q) = , 1 • (IH.3) 

At NLO, only the additional parameters CQ S / t enter. In Z-parameterization, they are chosen 
such that the residues of the poles in the dibaryon propagators do not receive any corrections 



beyond NLO; see again Ref. 



31| for details: 



M 

cos/t = —„ {Z s /t ~ 1) , (IH.4 

where Z s u = 1/(1 — j s /tPs/t) is related to the effective range p s u. Up to NLO, the dibaryon 
propagators in Eq. . 3h are modified by an insertion of the effective-range term c 0s / t of 
the Lagrangian of Eq. flHI.lj) : 

^ +NL °(<7o, q) = D™(q , q) + D^°(q , q) , (III.5) 
with the NLO correction (Z-parameterization variant in the second line) 



D^°(qo, q) = iD$(Q>, q) (-<*./*) ( * " ^ + ^ ) *°#fa>> q) (III. 6) 



%/t + \/^-Mq -ie Z s/t -l 



ls/t-^ g T -Mq -ie 



2%/ 



Calculations with a deuteron as an external state require wave function renormalization, 
Z t . In Z-parameterization, the LO expression is 



Z\° = ^ , (IIL7) 



and up to NLO 



z LO+NLO = 27* . 2jt 1 ( } 

* M M (l-j t p t ) 

With Zt = 1.6908, the NLO correction results in a 70%-shift from the LO value. While 

this contribution is much larger than expected from naive power counting, there are no 

further corrections to Z t at higher orders. This has the important advantage that the 

correct asymptotic normalization of the deuteron wave function at large distances r is exactly 

reproduced at NLO, with no corrections at higher orders: 



^dcutc r on(r oc) = . (III.9) 

V Zti r 

Taking into account the unusually large NLO term therefore significantly increases overall 
convergence of the expansion of EFT(^), as demonstrated e.g. in Refs. 



31 



33]. 



7 



B. Three-Nucleon Sector 



The consistency requirement of including a three-nucleon interaction (3NI) in the 2 Si 
channel even at LO is discussed in reviews; see e.g. Refs. 18l-l20|. The corresponding Lag- 
rangian is given by 



C 



3JV 



y 2 M Hq(A) 
3A 2 



[dl(a t N) - df(r A N)] j [dHa.N] - df(r A N)} 



(111.10) 



where Hq(A) denotes the three-nucleon coupling, which depends on the UV regulator A. This 
is the only parameter of Nd scattering up to NLO not determined from NN experiments. 
The 3NI strength Hq(A) can be chosen to reproduce the triton binding energy, or the 2 Si- 
scattering length. Choosing different low-energy data to fix it provides one method to 
estimate the theoretical uncertainties in Sec. IVI El 

The parity-conserving nd scattering amplitude is found by solving a Faddeev equation, 
see e.g. Ref. 31] . Its pictorial representation in Fig. [T] specifies the center-of-mass kinemat- 
ics: The total non-relativistic energy is E, and the momentum of the incoming (outgoing) 
deuteron is k (p). The on-shell point is at E = |£ — jf + ie and p = k. 





(2M' ^) (ILp p) 



Figure 1: Nucleon-deuteron scattering at LO. Shaded ellipse: LO amplitude t LO ; thick line: LO 
dibaryon propagator V LO of Eq. (jIII.12p : thin line (/C): propagator of the exchanged nucleon; disc: 
PC 3NI (H%°). 



In the dibaryon framework, the three-nucleon system receives contributions from Nd s and 
Ndt configurations, which are conveniently taken into account by decomposing all operators 



and amplitudes into the so-called cluster-decomposition basis (see Ref. 
our App. IA1 . 
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, App. A.l] and 



O = N, 



t 



( 



0(Nd t Nd t )\ 0(Nd s -)• Nd t ) j A 



b6 



\0{Ndt ^Nd s )f 0(Nd 



Nd s ) B A 



4 



N a 



(III.ll) 



s 



where the spinor and isospinor indices (a,/3) and (a, b), respectively, are often suppressed 
in the following. The vector indices i,j = (1,2,3) label the spin Pauli matrices while 
A, B = (1, 2, 3) label the isospin Pauli matrices. As an example, the LO dibaryon propagator 
in the cluster-decomposition basis is defined by a diagonal matrix in terms of the dibaryon 
propagators Eq. ( 1III.3I) : 



/ 



V LO (q ,q) 



\ 



(111.12) 



Df°{qo,q) 
y D^°(q ,q)J 

An analogous expression holds for the NLO correction of Eq. ( 1III.6I) . 

The spin-quartet channel only receives contributions from Nd t — > Nd t , and the cor- 
responding amplitude t q is the solution to an integral equation that only involves the (11)- 
element of the cluster matrix (the argument (E;p iQ ,p out ) applies to each entry in the matrix): 

A 

(E; k,p) =-4irJC {L) (E;k jP ) 



f t? ( 



v° °J 



dqq 2 )C^(E;q,p) 





V LO (E;q) 


(al) 

hq 










V° V 









{E;k,q), 



(111.13) 

with A the UV regulator. The projection of the exchange-nucleon propagator onto a specific 
orbital angular momentum L is 

Pl(cosA) 



1 r 1 

K {L \E-q,p) = - dcostf 

2 J-i 



p 2 + q 2 — ME + pq cos 9 



-l) L fp 2 + q 2 -ME\ 
-Ql 



pq 



pq 



(111.14) 

with Pl{ z ) and Ql(z) the Lth Legendre polynomials of the first and second kind with 
complex argument, respectively |34|, and 9 = Z(p; q). 

In the spin-doublet channel, amplitudes with different cluster-decompositions mix, so 
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that with t^ly denoting the amplitude for the process Nd x —> Nd y and x, y G {s,t}: 



AL) AL) 
l dM L d,st 



AL) AL) 
\ L d,ts L d,ss J 



(E;k,p) = 2ix 



dq q 



^ L) (E;k,p) 



JC {L) (E;q,p) 



xV LO (E;q) 




2H (A) 
A 2 



^ a l2^o(A) 



'0 



A 2 




(E;k,q) . 



(111.15) 



Since there is no partial- wave mixing in the PC sector even at NLO, the quartet- and doublet 
scattering amplitudes can be combined into one cluster matrix 

A 



/ ti L) 



( AL) AL) 
L d,tt L d,st 



(E;k,q) 
\ 



for the spin-quartet, S = § 



(111.16) 



AL) AL) 
L d±s b d,ss 



(E; k, q) for the spin-doublet, S 



where S is the spin, L the orbital angular momentum and J the total angular momentum 
of the 2S+l Lj partial wave considered. 

For the NLO PC amplitudes, we use the so-called "partially-resummed" formalism, in 
which the kernel and inhomogeneous part of the integral equations are expanded to NLO 
and then iterated, see Fig. [2] and Ref. 35(. This modifies the dibaryon propagators in 
Eqs. 13[4TlI. 15[) by an insertion of the effective-range contribution, i.e. by replacing 
V LO (E;q) with x> LO+NLO (E; q), see Eq. fllTLoT) . No new PC 3NI enters at NLO, but in 
order to reproduce the three-nucleon observable, Hq has to be adjusted at NLO. This is 
most conveniently achieved by dividing H into a LO piece Hq and a term at NLO, H^ . 
In this approach, in addition to all LO and NLO contributions, some higher-order contribu- 
tions are also included in the amplitude referred to as t LO + NL °. This does not increase the 
accuracy of the calculations, which is still set by the order to which the kernel is expanded. 
Figure [3] shows the cutoff dependence of the 3NI H at LO and LO+NLO in the partially- 
resummed formulation. While Hq varies considerably, observables are cutoff- independent, 
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Figure 2: Top row: Nucleon-deuteron scattering in the partially-resummed formalism at NLO. 
Middle row: all LO and NLO contributions. Bottom row: some of the diagrams which are strictly 
higher order than NLO but resummed for convenience. Cross-hatched ellipse: t LO + NLO as defined 
in the text; cross: insertion of the effective-range correction Co s /j to the dibaryon propagator; 
diamond: PC 3NI H^ L °. 





j < 




♦ 


i 1 


\ - 
f \ : 
\ ♦ • 



200 500 1000 2000 5000 

A [MeV] 



Figure 3: Cutoff dependence of the 2 Si-wave PC 3NI H (A), determined to reproduce the 2 Si 

2 2 

scattering length. ♦ (dashed lines): LO; • (solid lines): LO+NLO in the partially-resummed 
formulation. The linear extrapolations are only meant to guide the eye. 



see e.g. 3l| . This point will be of importance in the discussion of renormalization of the PV 
amplitudes in Sees. IVIB[ IVI CI and App. [Bj 
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The following parameters are used 3l|] : he = 197.327 MeV fm, the isospin-averaged 
nucleon mass M = 938.918 MeV; j t = 45.7025 MeV, 7, = -7.890 MeV, Z t = 1.6908, and 
Z s = 0.9015 are the effective-range parameters of the NN system; and the nd 2 S 1 scattering 
length 03 = 0.65 fm j^j] or triton binding energy B 3 = 8.48 MeV determines the 3NI H (A). 



IV. PARITY- VIOLATING LAGRANGIAN 



The leading-order PV Lagrangian relevant to our calculation is given in terms of five S— P 
wave transitions 



C 



PV 



+ 9 < lAiZ 3 o) 0) df (N T a 2 a-r 2 r A idN) 
+ 9(AiZl) 0) df (N T a 2 a ■ r 2 r B dN 
+ 2( ( 2: 2 T 0) Z AB df (N T a 2 a ■ r 2 r B Wn) 
+/ 5l " 3pi ) e ijk d? (N T a 2 a k r 2 T 3 d^N 



H.c. + . . . 



(IV. 1) 



where aOdb = aOdb — (da)Ob, O is a spin-isospin operator, and X = diag(l,l,— 2) 
is a dia gon al matrix in isospin space. For equivalent Lagrangians in different bases see 
Refs. 



diagonal 



We note as an aside that the relation 
non-dibaryon couplings C( x_y ) of Ref. 



Detween these PV dibaryon couplings g( x Y ^ and the 
l| differs from that given in Ref. p] due to different 



conventions used in the PC sector. The general expression remains 
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(X-Y) 



A 



(IV.2) 



with A s for X = x So and At for X = 3 Si, but the values for A s u and y in the Z- 
parameterization used here differ from the conventions of Ref. 2]. 

Higher-order contributions to the Lagrangian of Eq. (lIV.ip are suppressed by additional 
powers of Q. Corrections to S— P wave operators are expected to be suppressed by Q 2 
since they contain the same spin-space structure as the terms in Eq. ( lIV.ip . but with two 
additional derivatives. The effects of different partial waves, such as P— D wave mixing, are 
suppressed even further. 

As in the parity-conserving case, a simplistic application of the power counting suggests 
that parity-violating 3NIs first start to appear at N 2 LO. Unlike the PC case, this simplistic 
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power counting is valid for PV nd scattering; parity-violating corrections to S— P wave 
transitions from P V 3NIs do not contribute at LO or NLO . The Lagrangian of Eq. fllV.lj) 
is therefore sufficient to determine PV nd scattering up to and including NLO. 



V. NEUTRON-PROTON SPIN ROTATION 



In the dibaryon formalism, the only diagrams contributing to np spin rotation at LO 
are the tree-level diagrams shown in Fig. HJ The corresponding non-zero amplitudes are 




Figure 4: LO PV two-nucleon diagrams. Square: parity-violating two-nucleon vertices. 



obtained by using the LO dibaryon propagator of Eq. (1III.3I) : 

( 3 £i— 3 Pi) 

LW[nt p|->- nt pt] =i8V2ky - — , 

It + ik 

9 ( 3 Si- 3 Pi) 

iM [ni pi-)- n| p|l = - i8V2ky — , 

It + 1« 



iM[v\ p|^ n| pi] =iAV2ky \ + 
\M [ni pt^ ni pt] = - \A\Fik 



7t + ls + ik 

9 [bl Fl) , 9 {A I=0) - Z 9 (AI=2 ) 



It +ik % + ik 

with k the magnitude of the center-of-mass momentum. 

Applying Eq. ( 111.81) with \i = M/2 for the reduced mass and /ci a b = 2k, the spin rotation 
angle for polarized neutrons on a hydrogen target is given by 

^^f = i^ E Re[M + K)-M_K)], (V.2) 

m p =± j 

where m p = ±~ is the (initial and final) proton polarization. Since the (thermal) external 
nucleon momentum k appears only in the dibaryon propagators, it is negligible compared 
to the parameters j s / t . With Eq. ( IV. lj) . the PV rotation angle for a hydrogen target at LO 
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is given by 



i 



p dZ 



--2V2nM 



29 



(Si Pi) +g (s 1 # (A7=0) -2g (AI=2) 

Is 



LO 



1/ 



3.4 ( 2/*- 3 ^ + - 19.5 ( - 2^ 



(V.3) 



rad MeV"2. 



where the g( x Y ^ carry units of MeV 3 ' 2 . 



/ CQs/t \ 




Figure 5: NLO PV two-nucleon diagrams. 



Since no new PV operators enter at NLO, as discussed in Sec. IPV[ the scattering amp- 
litudes at NLO are given by the diagrams in Fig. [5j Using the dibaryon propagator of 



Eq. (lIIL5j) . the LO+NLO result is obtained by multiplying Eq. fTVlj) by — ^ — for each 



term - x +ik if the neutron momentum is neglected compared to j s /t. This correction is 
35% for the P— 3 Si waves and 5% for the P— x So waves, in agreement with the naive power 
counting estimate of the EFT(^) expansion. Our expression up to NLO becomes 



1 



dl 



LO+NLO 



l[4.5±0.5] (2 9 < 3s '- Sp '>+/ s '-' F '>' 

- [i8.5 ± i.9] (silit"' - ^S:? 1 ') ) rad MeV ~* 

where we have conservatively assigned errors of 0(Q 2 ) ~ 0.1. 



(V.4) 



VI. NEUTRON-DEUTERON SPIN ROTATION 

A. nd partial-wave amplitudes at leading order 

The parity-violating nd scattering amplitude receives contributions from tree-level 
(Fig. E]), "one-loop" (Fig. [7j), and "two- loop" diagrams (Fig. EJ). This nomenclature refers 
to the number of loops that contain a parity-violating interaction; the strong amplitudes 
receive contributions from an infinite series of multi-loop diagrams. All graphs use the same 
interaction kernel, namely the tree-level PV diagrams in the off-shell kinematics specified 
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(a) 



Figure 6: LO tree-level PV diagrams and off-shell kinematics for convolution in Figs. [7] and [8l 





Figure 7: LO "one-loop" PV diagrams. "Time-reversed" contributions (diagrams obtained by 
reading above ones from right to left) not displayed. 




Figure 8: LO "two-loop" PV diagrams. "Time-reversed" (as defined in Fig. [7]) contributions not 
displayed. 



in Fig. |6j These are in turn computed using the Lagrangians of Eqs. (IHI.lj) and (1IV.1I) . 
and with the propagator /C of the exchanged nucleon in Eq. (IIIL14I) as well as the cluster- 
decomposition version of the dibaryon propagator in Eq. ( 1III.12I) . Amplitudes between 
two states with angular momenta L and L' scale as k L+L ' for small nucleon momenta, see 



e.g. 



33, p. 381 Eq. (157)]. Since spin-rotation measurements are performed with neutrons 



with momenta k 1 MeV, only transitions between S and P waves contribute at this order, 
i.e. L, V e {0; 1}. The projectors necessary to filter out specific three-nucleon partial waves 
are constructed and listed in App. [A] 
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While the PV Lagrangian contains five independent couplings, forward Nd scattering at 
low energies only depends on three independent (isospin-dependent) linear combinations: 

S 2 = 3/*-*> - r 3 /^ , (VI.l) 

1 - 6 9 {A I=0) + ZT 39(AI=1) ■ 

Since the Nd system is an iso-doublet, the PV coupling g^i=^ cann °t contribute. Here we 
are interested in scattering on a neutron, so only the (22) component of the isospin matrix 
is needed, effectively replacing r 3 with —1. 
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The projected scattering amplitudes of Fig. [6] (a) are: 
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(VI.2) 
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with E = ||^ — ^ + ie as before and 



+ g 2 - ME 
pg 



(VI.3) 



The structure of these matrices in the cluster space of Eq. (IHI.lip already appeared in 
Ref. [33]. The projector {Q r s ) a p = 5 r s 5p — \{(T r a s ) a ^ onto the spin quartet is defined in 



Eq. (IA.3I) . The index a (f3) denotes the spin of the outgoing (incoming) nucleon, and r (s) is 
the spin component of the outgoing (incoming) 3 Si dibaryon. The diagrams of Fig. Mjo) with 
the PV interactions on the lower line can be obtained from those of Fig. |6]^a) by Hermitian 
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conjugation in cluster, spin, and isospin space: 

M( b) lX^Y;q,pY a s p=(M {a) [Y^X;p,q] sl5 ra y , (VI.4) 

where X (Y) denotes the partial wave of the incoming (outgoing) state and the spin indices 
are made explicit. To find the contributions of Fig. [6] to PV scattering, choose the on-shell 
point p = q = k: 

LMtree[^ ^Y;k} = \M {a) [X -> Y] k, k] + iM {b) [X -> Y; k, k] (VI.5) 

The "one-loop" and "two-loop" contributions of Figs, and [8] are generated by convoluting 
the PV tree-level results with the PC amplitudes of Sec. IIIIl following the calculation of 
higher-order corrections in the PC sector in Ref. 29]. 

When the PC amplitude is attached to the left of the PV kernel in a "one-loop" diagram, 
as in Fig. [TJ the q integration picks the nucleon pole, q = — + ie, and the angular 
integration is trivial, 

A 

iMwoop, pc leftt* ->Y ] k]=Y J ^ [X -> Y; q, k] W LO (E; q) it LO [X; k, q] . 

(VI.6) 

Note that we choose the same UV regulator A as in the integral equations for the PC 
amplitudes, Eqs. ( 1HL13|/1TII.15I) . When the PC amplitude is attached to the right, the 
amplitudes are obtained by reading Fig. [7] as if the time direction were reversed: 

A 

iMi-?oo P , PC rights ^ Y;k] = J ^ ^°\y; K q] iV LO (E; q) iM® [X Y; k, q] , 

j=a,b o 

(VI.7) 

where we used that the PC amplitudes are time-reversal invariant, i.e. identical when ex- 
changing incoming and outgoing nucleon momenta, t[X;p,q] = t[X;q,p\. The "two-loop" 
convolutions of Fig. |8] are 

iM\Z op [X^Y ] k}= (VI.8) 

A A 

E/^ / '^^;M^;p)^[^ ^Y-q,p]iv^(E ]q )it^[x-k,q} . 

j=a,b 

We numerically solve the integral equations for the PC amplitudes t, Eqs. ( 1III.13| / UH.15I) . by 
the Hetherington-Schick method 40| in a Mathematica code, with A a hard cutoff corres- 



ponding to a step function, as detailed in Ref. 



311 ] . The convolutions for the PV amplitudes 
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Eqs. (1VI.6| / fVL7) / 1VI.8j) are performed numerically, using again a Mathematica code. The 



numerical error of < 0.1% of these procedures is negligible relative to the systematic uncer- 
tainties discussed in Sec. IVI El 

Nucleon-deuteron scattering corresponds to the (11) element of the cluster matrix. Mul- 



tiplication with the LO wave-function normalization factor a/ Z\° , Eq. (lIII.7p . for each ex- 
ternal deuteron leg results in the renormalized, physical scattering amplitude Ai R ° between 
partial waves X and Y at nucleon momentum k: 



iM L R °[X ^Y;k] = (Jz™, 0) (iM trcc [X ^Y;k] + iM™ op> PC right [X Y; k] 
+ i^i-?oo P , pc Mtl* ^Y-k} + iM\° loop [X -> Y] k]) ( 







r 3 ->-l 

(VI.9) 



B. nd partial- wave amplitudes at next-to-leading order 



No new PV operators contribute to the NLO amplitudes since higher-order PV interac- 
tions are suppressed by at least Q 2 as described in Sec. IIVI There are three types of NLO 
corrections. The first type is given by diagrams with one insertion of an effective-range 
correction to the dibaryon propagators proportional to c 0s / t , see Eq. (IIII.4I) . At this or- 
der, the momentum-independent PC 3NI parameter H Q (A) requires an additive adjustment, 
called Hq LO (A), to recapture the low-energy 3N observable that initially fixed Hq(A) at 
leading order. This results in a second type of NLO correction to the nd system: diagrams 
with H^ LO (A) inserted once between LO amplitudes. The third type of corrections comes 
from the change in the dibaryon wave-function renormalization and is taken into account 
by replacing Z\° in the LO PV amplitude of Eq. (TVL9]) by ^lo+n^ see Eq. f l!TL8|) . 



In a strictly perturbative approach, there are two classes of diagrams at NLO. The class- 
I diagrams of Fig. [9] are generated when the PC LO amplitudes t LO [X;k,q] are replaced 
with the partially-resummed amplitudes t LO+NLO [X; k, q] in the topologies of the LO PV 
"one-" and "two-loop" diagrams of Figs. [7] and El i.e. t LO [X;k,q] -> t LO+NLO [X; k, q] at 
each occurrence in Eqs. (jVI.6| / fVT77| / 1VI.8[) . see Fig. [TUJ In addition to the NLO class-I 



diagrams, the amplitudes obtained from Fig. [10] also contain all LO contributions as well 
as a partial resummation of higher-order effects. In analogy to the notation at LO, the 
resulting contributions are called A4i_?o^ L pc lefT^ ^ ^i-?otp L pc right I [ x ~> Y ', k }i 
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CQs/t PV Co s /t PV 




M LO PV # n NLO PV # n NLO PV tf n NLO PV 




H^ LO PV tfn NL ° PV M LO PV ff n NL ° PV 




Figure 10: Reformulation of the class-I diagrams of Fig. [9] using the amplitudes t LO + NLO of Fig. El 
This also replaces the LO PV contributions of Figs. [7]and[8j Diagrams with the PV vertex on the 
lower line as well as "time-reversed" contributions not displayed. 



and M 



LO+NLO, class-I 
2-loop 



[X — > Y- k]. This approach has computational advantages since it avoids 



diagrams with up to four numerical integrations and convolutions with parity-conserving 
LO full-off-shell amplitudes t LO \ 2S+1 L k p . q], p, q ^ k. In addition, the diagrams in class I 



contain divergences, studied by Refs. 



38 



41 



42[, that in a strictly perturbative calculation 



are cancelled by corresponding contributions in class-II diagrams, to be discussed shortly. 
These divergences are challenging to treat numerically, but are avoided in the approach used 
here. The amplitudes corresponding to the diagrams of Fig. [10] are individually renormalized; 
each diagram approaches a unique and finite limit as the cutoff is removed, A — > oo. 
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Diagrams that are not generated by the replacement t LO [X;k,q] — > t LO+NLO [X; k, q] 
in Figs. [7] and |8] are referred to as "class-II" diagrams and are shown in Fig. [TTJ The 



cos/tPV co s/t PV PV 




CQs/t 

Figure 11: Class-II diagrams. As discussed in the text, the PC amplitudes are taken as £ LO + NLO 
for 2 Si and 4 S3, and as t LO for 2 Pj and 4 Pj. Diagrams with the PV vertex on the lower line as 



2 2 

well as "time-reversed" contributions not displayed. 



corresponding amplitudes are given by 

A 

^i° ^f cit [X -> Y ; *] = £ | ^£ Wl® [X -+ Y- q, k] W^°(E; q) it[X; k, q] 



(VT.10) 

iM^^ljX ^Y;k]=J2 [ m k, q] ^°(E; q) \M® [X -+ Y; k, q] , 



A 

,2 



]=a,b o 



(VI.11) 

A A 

E/|^ J^mk,p}^°(E;p) i/A^ [X — >■ Y; q,p] iD NLO (E; q) 

j=a,b 

+ iV NLO {E;p) iM {j) [X -»• Y;q,p}iV LO {E;q)}it[X;k,q} . 



Since the kernel of these convolutions is now already NLO, it appears that the LO PC 
amplitudes t LO [X; k,q] can be used directly. This is the strategy implemented for the PC 
P-wave amplitudes. 

For a convolution of the class-II diagrams involving the LO PC 2 Si-wave amplitude 
the situation is more complicated. It was demonstrated on general grounds in Ref. jssj], 
and confirmed numerically here that when using t h ° the "one-loop" amplitudes diverge 
logarithmically, and the "two-loop" amplitudes diverge as g - 23 --; see App. |B] These are 
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In a 



not divergences that are removed by additional parity-violating 3NIs at NLO 
strictly perturbative calculation in the PC sector, the 3NI Hj? LO absorbs a linear divergence 
generated by insertions of the effective-range term 4l|, |42|. It is this linear divergence 
in Hq LO which, when inserted next to a PV interaction (see class-I diagrams in Fig. |9]), 
would renormalize the divergence in the class-II contributions. However, in the partially- 
resummed formalism used here for the class-I diagrams //LO+nlo ^ oes no ^ diverge linearly 
as A — > oo, see Fig. [3J since the high-off- shell momentum part of the scattering equation, 
and hence of the amplitude, becomes softer. Therefore, the treatment of class-I diagrams 
using the partial resummation technique, see Fig. \TU\ removes those divergences from the 
class-I diagrams. This in turn means that the class-I diagrams do not renormalize the 
divergent class-II contributions, so that the divergences in the class-II diagrams require 
separate renormalization. By using the renormalized (partially-resummed) PC amplitude 
t LO+NLO [ 2 Si; k, q] for the 2 Si channel the UV behavior of the class-II diagrams is changed 
and no divergences appear. As demonstrated below in Fig. [T2l this leads to renormalized, 
cutoff-independent PV amplitudes. We therefore choose the renormalized PC amplitude 
t LO+NLO [ 2 Si ; k, q] for this channel. 

We also choose the LO+NLO version for the PC 4 S3-wave amplitude in the class- 

2 

II diagrams. This is not required by renormalization, as the loop-integrations in 



Eqs. ( 1VI. 10| / fVI. 1 l| / fVT. 1 2j) converge well in this channel, see again Ref. 38j. However, neutron 
spin-rotation in deuterium is particularly sensitive to the scattering lengths of the nd system 
since it is a process essentially at zero kinetic energy. It is therefore important to reproduce 
the experimental values of the scattering lengths. In particular, the 4 S3 scattering length is 
a factor of 10 larger than the 2 Si one, so it will most likely dominate spin-rotation observ- 
ables. At LO in the Z-parameterization used here, the 4 S3 scattering length of 5.1 fm would 
require a 25 percent correction to achieve the experimental value of [6.35 ± 0.02] fm J36], 
while the NLO amplitude leads to 6.4 fm, in close agreement with the experimental value, 
see Ref. [3JJ. We therefore use the NLO PC amplitude t LO+NLO [ 4 S3 ; k, q] for this channel. 



This choice does not violate the power counting of the theory and yet improves the utility 
of our result. 
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To summarize, the PC amplitudes used in Eqs. (1VI. 10[ / TVI. 1 1|4VT. 12|) are: 



f f LO+NLO k g l for X e {2g 4 Ss l 

t[X;k,q} = { 5 ' 5 ' (VI.13) 

t LO [X;fc,g] for X G { 2 P; 4 P} . 



i 



Further details and comparisons to other choices are discussed in App. [B] 

The renormalized LO+NLO PV amplitudes are obtained by summing and multiplying 
with the LO+NLO wave-function renormalization: 



^LO+NLO 



[X ->Y;k]= (y^ LO+NLO , o) (iM tTCC [X -> Y; k] (VI. 14) 

+ iM™£&'£r I [X ^Y-k]+ iM^ v %l h [X -+ Y; k] 
+ ^SS'II —^Y;k}+ iM^^lJX -+ Y; k] 



LO+NLO 



, - a /LO+NLO, class-Ir v > v n , ■ a y*NLO. class-IIr v > v ;Af v^t 

+ 1 - M 2-ioo P t x ->• Y ~A + ^2-ioop [*-> y ;fc]H 

Finally, we reiterate that these amplitudes are complete up to NLO, but also contain some 
higher-order contributions. The PC amplitudes £ LO+NLO used are partially resummed. In 
addition, we choose to multiply the entire LO+NLO amplitude by the LO+NLO wave- 
function renormalization, rather than the LO PV amplitude -M^ by ^ LO + NLO an d the 
NLO correction (.Mp 0+NL0 — .Mr ) only by Z L °. It has been demonstrated repeatedly that 
this speeds up convergence since the larger-than-usual wave-function renormalization cor- 



rection at NLO is expected to be t 
operators; see e.g. Refs. 
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31 



re dominant correction to insertions of higher-dimension 



33]. 



The computational effort can be halved by taking advantage of the fact that the individual 
tree-level amplitudes Ai^ and Ai^ are Hermitian conjugates, Eq. ( 1VI.4I) . and that the PC 
amplitudes are time-reversal invariant, t[X;p,q] = t[X;q,p\. The amplitudes of partial 
waves X — >• Y and Y — > X are therefore related: 

M l-loop, PC right 

[X ^Y;k} = (M l-loop, PC left [Y^X-k}) ] , (VI.15) 

Afi-ioo P [V ^Y;k] = (M 1Aoop [Y -> X- k}) ] , (VI.16) 

-M 2 -ioo P [V ^Y-k} = {M 2 _ Xoop [Y -+ X- k}) ] , (VI.17) 

M R [X -> Y; k] = M R [Y -> X- k] . (VI.18) 
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C. Cutoff independence of partial- wave amplitudes 

Numerically, it is difficult to perform some of the integrations at very low energies. How- 
ever, between k = 0.01 MeV and 10 MeV the amplitudes of the nd system deviate from the 
linear relation only at the percent-level. For the spin rotation angle, the relevant quantity 
is M R /k, see Eq. (fIT8]) . The difference in M R /k at k = 1 MeV and 0.01 MeV is less than 
0.3%. All results quoted are for k = 1 MeV. 

For the analysis of the cutoff dependence of the partial-wave amplitudes, we decompose 
them in terms of the PV interactions <S»i, S2, and T of Eq. IVI. II 

k 

+ d[X -> Y; 5 2 ](A) S 2 + d[X -> V; T](A) T . 

In Figs. [12] and [131, the dependence of the functions d on the cutoff A used in the PC 3N 
integral equations Eqs. (1III. 13)4111. 15j) and in the convolutions of Sees. IVI Al and IVI Bl is 



shown for each partial-wave and order. Since a unique, finite limit exists as A — > 00, the 
result is properl y re normalized in each partial wave. These results confirm the theoretical 



findings of Ref. |38[ that no PV 3NI is required at LO and NLO. Comparison with Fig. [3] 
also shows that at both LO and NLO the physical amplitudes are smooth even where the 
3NI H diverges. 

A small kink in the cutoff dependence of the partial wave amplitudes is seen at NLO for 
A 3000 MeV. It is evident in all partial waves, and so is not related to the renormalization 
of the 3NI H of the 2 Si-wave. The phenomenon appears only at cutoffs well beyond the 
breakdown scale of 200 MeV and hence does not have any impact on our final result. 



D. Translating partial-wave amplitudes into neutron-spin rotation predictions 

For computational convenience, the calculations so far have been performed in a partial- 
wave basis. However, in order to obtain the spin rotation angle, we need to determine 
the forward scattering amplitude for specific neutron helicity states. Recall that in the 
kinematics of Sec. IHIl the incident and outgoing nucleons have momentum —k = —ke z , so 
that an incoming neutron N aa (—k) with positive helicity corresponds to choosing the spin 
and isospin components a = 2, a = 2. The relation to the results in the partial-wave basis 
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Figure 12: Cutoff dependence of the functions <i[ 2 Si — > Y; coupling] (A), Eq. (|VI.19p . Columns 

2 

(left to right): partial waves 2 Si — > 2 Pi, 2 Si — > 4 Pi; rows (top to bottom): coefficients of 5 l5 

2 22 2 

52, T. * (dashed lines): LO, 3NI Hq determined from nd 2 Si-scatt. length; A (dashed lines): 

2 

LO, 3NI from triton binding energy; • (solid lines): NLO, 3NI from 2 Si-scatt. length; ■ (dotted 

2 

lines): NLO, 3NI from triton binding energy. The linear extrapolations are only meant to guide 
the eye. Notice the different scales on the vertical axes. Amplitudes which are identically zero are 
not displayed. 
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Figure 13: Cutoff dependence of the functions <i[ 4 S3 — > Y; coupling] (A), Eq. (IVI.19I) . Columns 

2 

(left to right): partial waves 4 S3 — > 2 P3, 4 S3 — > Pi; rows (top to bottom): coefficients of Si, S2, 

2 22 2 

T- * (dashed lines): LO; • (solid lines): NLO. Comments as in Fig. [T2l 



is obtained by inserting a complete set of projection operators: 

(d,N bp (-k)\M\d l N aa (-k)) = J2(d,N b ^-k)\Y)(Y\M\X)(X\d l N aa (-k)) , (VI.20) 

XY 

where (F|A^|X) = .Mr.[X — >■ Y] are the partial-wave amplitudes calculated above. The 
sum over X, Y stands for all partial waves including spin and isospin polarizations. 
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Using the three-body projectors V[X] constructed in App.|A]from Eq. (1A.14|) . the partial- 
wave projected matrix elements of the neutron-deuteron state are 

(X = { 2S+l LrM{M}) \d l N aa (-k)) = [{V^ M \[ 2S+1 Lj}f aa \ n , (VI.21) 

with appropriate vector and spinor indices {M}. The outer bracket [.]n indicates that only 
the (ll)-entry of the cluster-decomposition matrix is physically allowed. As an example, the 
projection onto the 4 Pi-wave is 

'3 



'PildiN^i-kek)) 



3 \0 



a u a i 



(VI.22) 



where the indices 0, b are needed to specify the spinor and isospinor magnetic quantum 
numbers of the 4 Pi state, and = (0,0,1) for forward-scattering. The positive-helicity 
neutron-deuteron amplitude is then: 



M + = Y,M K [X->Y] 1 - £ ((-P{M} iA [Y]y V 

XY ielO;±l> L 



\<r>{M} 



i A 



[X] 



/3=2,a=2;a=2,b=2 



. (VI.23) 



n 



iS{0;±l} 

Note that the iso-vector index A is irrelevant for the final (ll)-component of the cluster- 
decomposition matrix. This formula can also be understood as a decomposition into the 
reduced matrix elements, represented by the amplitudes A^r[A — » Y] calculated above, 
and the equivalent of Clebsch-Gordan coefficients weighting the various partial waves. The 
negative-helicity amplitude is M._ = Using the partial-wave projectors of App. |A] 

leads to our final result for the neutron spin-rotation in deuterium from Eq. flll.8j) with 
fi = 2M/3 for the reduced mass: 
1 2M 4 



p dl 



6k 9 



Re 



M R [ 2 Si ^ 2 Pi;k] -2v / 2A^ R [ 2 Si ^ 4 Pi;A;] 

2 2 J L 2 2 ' 

-4M R [ 4 S 3 2 P 3 ;/c] -2v / 5A4 R [ 4 S 3 -> 4 P 3 ;£;] 



(VI.24) 



where the amplitudes for P to S wave transitions are taken into account by using the relations 
of Eq. (jVEEgp . 



E. Numerical nd spin rotation result and error estimates 

For the detailed discussion of theoretical errors below, we decompose the spin-rotation 
result of Eq. (I VI. 24 j) in terms of functions c[(X — Y)](A) which multiply the PV interactions: 

~ ^ = c[(% - x Pi)](A) / Sl " 1Pl) + r 3 c[(% - 3 Pi)](A) / 5l " 3pi) + c[7l(A) T (VI.25) 
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The isospin matrix T3 is replaced by —1 for neutron spin-rotation, and T = ^9^°Z^ + 
2 r 35 , (Ai=if ' ) ' see Eq- (fVI-iP - Figure [TJ] shows the cutoff dependence of the functions 
c[(X — Y)](A). As expected from the discussion of the cutoff dependence of the partial- 
wave amplitudes, the result for the spin-rotation angle is also properly renormalized. 
The final result for neutron spin-rotation in deuterium is: 

I = (j 8 ± 8] g es^ Pl ) + [170 ± L7] g e Sl -^) (VL26) 

+ [2.3 ± 0.5] (3^2?-^ - 2g^) )rad MeV- , 

where the estimate of the residual theoretical uncertainties is justified below. It will turn 
out that this estimate is rather conservative; drawing from the experience in calculations 
of PC 3N observables using the Z-parameterization, theoretical systematic errors might be 



estimated to be on the order of 3%, see e.g. 3l|, |43[. For now, we note that the theoretical 



errors of 10% to 20% are comparable to the statistical and systematic errors expected of the 
most ambitious planned experiments. 

We use three methods to estimate theoretical uncertainties, with the EFT(^) parameter 
Q ~ |_^ s a conservative value for typical momenta in the Nd system on the order of j t (see 



Ref. (ai|). 



(i) At NLO, higher orders should contribute corrections of order Q 2 m 0.1 to the total 
result. 

(ii) Since EFT(^) is not valid at high momentum, a residual dependence of observables on 
momentum modes above the breakdown scale A ~ 200 MeV is an indicator of the size 
of higher-order corrections. Varying the momentum cutoff A in the Faddeev equations 
and convolutions of the 3N system from the breakdown scale to higher values provides 
an estimate of N 2 LO effects. The parameters g^ x ~ Y ^ of the PV Lagrangian Eq. (1IV.1I) 
are independent of the regularization scale A. Therefore, the c[(X — Y)](A) of the 
spin-rotation results of Eq. flVI.25j) must also be cutoff-independent up to higher- 
order effects. In particular, this implies that these coefficients should approach a 
unique, finite limit as the cutoff is removed. As in a variety of previous calculations 
in the 3N system, we consider cutoff-variations from A = 200 MeV to 5000 MeV, see 



e.g. Refs. |3jJ, 135 ]. 



2N 




Figure 14: Cutoff dependence of the functions c[(X — Y)](A) for the neutron spin-rotation angle 
in deuterium, (jVL25]). Notation as in Fig. [T2J Grey band: estimated theoretical uncertainties as 
described in Eq. (jVL26]). The linear extrapolations are only meant to guide the eye. Notice the 
different scales on the vertical axes. 



(iii) Finally, the size of higher-order effects can be assessed by using different low-energy 
data to determine the PC parameters of EFT(^). Since we chose Z-parameterization 
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precisely because of its well-established improved convergence, we do not vary para- 
meters of the NN system, e.g. by replacing (Z a u — 1) at NLO by the effective ranges. 
On the other hand, the strength of the PC 3NI Hq can for example be determined 
from the nd scattering length of the 2 Si-wave, or from the triton binding energy. The 
difference between both approaches is again a measure of N 2 LO effects, i.e. expected 
to be on the order of Q 2 ~ 0.1. 

Note that methods (ii) and (hi) do not apply to the 2N system. We will base our error- 
estimate for the nd spin-rotation coefficients on the most conservative of the above methods. 

The NLO corrections are as large as 70% relative to the LO result. The size of this 
correction stems from the unnaturally large residue of the deuteron pole, Z t — 1 ^0.7, in 
Z-parameterization. However, once this effect is taken into account, convergence at N 2 LO 
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33]. Different inputs to determine the PC 3NI affect 



and higher is actually improved 
only the 2 Si — > X partial-wave amplitudes. At NLO, this can lead to a change of up to 15% 
in the amplitudes -Mr[ 2 Si — > Y; k], in line with the error-estimate criterion (hi). 

In the functions c[(X — y)](A) multiplying the PV couplings g^ Sl ~ 1Pl ^ and g^ Sl ~ 3Pl \ 
Fig. [HI contributions involving the 2 S i-wave are doubly suppressed relative to those which 
contain the 4 S3-wave. Not only are their relative weights in Eq. ( 1VI.24j) small, they are 
also small in absolute size, as seen from the fact that the 2 Si-wave scattering length itself 
is a factor of ~ 10 smaller than the 4 S3 one. Therefore, applying criterion (hi) leads to 
very small variations. Varying the cutoff A G [200; 5000] MeV, criterion (ii), produces 
a range at NLO of c[( 3 Si - x Pi)](A) = [7.8... 8.1] rad MeV"^ and c[( 3 Si - 3 Pi)](A) = 
— [16.5 . . . 17.4] rad MeV~2. For error-estimate criterion (i), i.e. Q 2 ~ 0.1 of the total result, 
corrections are ±0.8 rad MeV" and ±1.7 rad MeV~2 ; significantly larger than the estimate 
from varying A. We therefore adopt the range from criterion (i) as a conservative estimate 
of the theoretical uncertainties in these coefficients. Overall, the amplitude .Mr[ 4 S3 — > 4 Pi] 
dominates c[( 3 Si — 3 Pi)](A), providing more than 80% of its total value. 

Due to symmetries, the function c[7~](A) does not receive any contributions from the 
4 S3 — > 4 P3 channel, and the amplitude .Mr[ 4 S3 — > 2 Pa] is very small. The 2 Si — > 4 Pi 
amplitude thus dominates c[T](A). Since this amplitude does not depend significantly on 
the input used to determine the PC 3NI, criterion (hi) significantly underestimates the 
theoretical uncertainty of the NLO calculation. Cutoff- variation, criterion (ii), is significant 
in this channel, mapping out the range c[T](A) = [1.8 . . . 2.8] rad MeV~2 ; while criterion (i) 
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provides an error-estimate of only ±0.3 rad MeV 2 . We therefore adopt the variation from 
the cutoff chana; conservative estimate of the theoretical uncertainties. While this 

error is ~ 20%, the magnitude of c[7~](A) is at most a third of the other two c[(X — Y)](A). 
If there are no accidental cancellations or significant differences in the magnitudes of the 
PV couplings themselves, the overall contribution of c[T](A) to the spin rotation angle is 
expected to be small. 



VII. ROTATION ANGLE ESTIMATES AND COMPARISONS 

Our results for np and nd spin rotation are given in terms of the LECs g( x ~ Y > of the 
Lagrangian flIV.lj) . As discussed in the Introduction, these couplings are not predicted by 



the EFT, but can be estimated on dimensional grounds. Since on the microscopic level 
the dominant PV process stems from a single weak gauge boson exchange, the PV NN 



interactions are suppressed relative to the PC ones by a factor ~ 7$r, where 74 is taken as 



M 



w 



a typical low-momentum scale. As in the PC interactions between two nucleons and the 
auxiliary dibaryon field in Eq. (lIII.2j) . an overall factor of \/ff should be included. In order 
to match the dimensionality [MeV -5 ] of the couplings g^ x ~ Y \ an additional inverse mass 
dimension enters. Since the PV couplings used here are renormalization-group invariant 
and encode short-distance physics, the mass scale can only be set by the breakdown scale 
A m w 140 MeV of EFT(^). Overall, therefore, we expect the magnitude of the PV 
couplings to be of order 

C*-y)| „ J** j ( Jl) 2 « 10-10 MeV-l • (VII.l) 



' M A \M W/ 

An estimate which is compatible with this number may be obtained by appealing to Ref. [l|, 
where a value for the combination of couplings involved in PV pp scattering at 13.6 MeV 



44l | was extracted. (Asymmetries measured at higher energies are not within the realm of 



validity of EFT(/).) Using this value as an estimate of all LECs, 

\g {x ~ Y) \ « 2 x 10 -11 MeV-i (VII.2) 



in our conventions. Note that these are only rough estimates. 

With these values for the magnitude of the PV couplings and a target density of p w 
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10 23 cm 3 , the magnitude of the spin-rotation signal in hydrogen of Eq. ( IV. 4 P is 

rad 
m 



,np 
^PV 



dl 



[10^ 7 • ■ • 10~ 6 ] ^ . (VII.3) 
Analogously, the magnitude of the spin-rotation signal in deuterium of Eq. ( 1VI.26I) is roughly 

[10-7... 10~ 6 ] — . (VII.4) 



dl 



We stress again that - without reliable values of the PV parameters - these results are 
dimensional, order-of-magnitude estimates and may well be off by factors of 10 or more. 

Another order-of-magnitude estimate can be obtained by using the parameter set of 
Ref. 2^1 extracted from the DDH "best estimates"; using these values (Eq. (31) of Ref. 26]) 
and adjusting for different sign conventions yields spin rotation values of 

^ rad 



iJip 



dl 



for np spin rotation and 



knd 



3 x 10" 



5 x 10 



in 



_ 7 rad 



(VII.5) 



(VII.6) 



dl m 

in the nd case. Note that there is a large spread in the DDH "reasonable ranges" surrounding 
these DDH "best estimates." 

We now compare our estimated result with the results from other calculations. Without 
measured PV parameters, none of these can be considered as more than order-of-magnitude 
estimates. Hence we do not normalize with respect to less than order-of-magnitude changes 
in choice of target density, for example. 

Our estimate for the np spin rotation angle ag rees well with earlier results obtained 



bj IM 



24j, which found rotation angles between 



in the DDH and hybrid formalisms [l^, 
5.15 x lO -7 ^ and 1.36 x 10 -6 ^, depending on model choices made for the couplings and 
the strong interaction potentials. 

The rotation angle in neutron deuteron spin rotation was determined in Ref. 15] using 
the DDH model and the Argonne vis interaction, supplemented with the Urbana-IX three- 
nucleon potential. The result is 



d(j> 



nd 
PV 



9.32 x 10" 



,rad 



dl m 
while a subsequent calculation also employing the At> 18 +UIX potentials found [1 



nd 
PV 



dl 



7.68 x 10" 



,rad 



in 



(VII.7) 



(VII.8) 
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This agrees with our findings that the coefficients of the PV couplings in np and nd rotations 
are of the same order of magnitude, and also agrees with our rough estimate of Eq. (jVTU]) . 
Note that result (I VII. 81) is just one of several nd spin rotation predictions provided by 
Ref. jl7| . who consider two other parameter set estimations in addition to a set with the 
"best values" from the wide ranges provided by DDH. These two parameter sets, collected 
by Bowman [45(, yield values of -6.82 x 10~ 7 and -8.91 x 10~ 7 — . 



In addition, Refs. 15| and [l7| performed a hybrid calculation in which the phenomen- 
ological wave functions of the PC sector are combined with a PV potential derived from 
EFT(^f). In this approach the consistency of treating interactions, wave functions, and cur- 
rents within a unified framework is lost. In particular, different degrees of freedom appear in 
the PC and PV interactions. In both references, results for two different values of the regu- 
larization parameter /i used in the PV EFT(fi) potential are given, while in the PC sector a 
particular pionful potential with a fixed parameter set is used. For the choices fi = 138 MeV 
and n = 1 GeV, the results for the coefficients multiplying the PV parameters differ by up to 
two orders of magnitude. While this regularization dependence can in principle be removed 
by running the PV couplings to absorb cutoff dependence, these large differences are an 
indication of the resolution mismatch between EFT(f) and the phenomenological poten- 
tials. Also note that different regularization schemes and degrees of freedom are employed 
for PC and PV interactions. Since the renormalization-scale dependence of the PV and PC 
couplings in the hybrid formalism is unknown, one cannot determine whether the calcula- 
tion is indeed cutoff-independent, and varying the cutoff cannot be used to assess residual 
theoretical uncertainties. On the other hand, as shown above, our results are independent 
of any choices related to regularization and renormalization within the errors estimated in 
the EFT approach. 



VIII. CONCLUSION AND OUTLOOK 



At NLO in parity- violating EFT(^), two- and three-body low-energy PV observables 
depend upon five unknown LECs, the g {x ~ Y) of Eq. flTV~ll . PV EFT(?f) is based on the 
symmetries of QCD and the weak interaction, and the power counting of the theory. In 
order to verify that QCD is appropriately encoded into this EFT, it must be demonstrated 
that the LECs extracted from independent observables are consistent and are of "natural" 
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size - 0(1) in dimensionless units. To overconstrain the system requires more than five 
model-independent calculations involving linearly independent combinations of the LECs, 
and the availability of the corresponding measurements. 

This paper provides two of these calculations: Using EFT(^) consistently in both PV 
and PC sectors of the nucleon interactions, we obtain model-independent results for neutron 
spin rotation in hydrogen and deuterium targets. At NLO they are given by 



1 d0 



PV 



p dl 



([4.5 ±0.5] (2/ ft " 3pi ) + /*-^)) 

-[i8.5 ± i.9] (« ) -2,fS::r° ) : 



rad MeV"2 



for np, and 



Lnd 



P dl 



([8.0±0.8] ^C 3 ^- 1 ^) + [17.0 ± 1.7] /*- 3 ft) 



[2.3 ±0.5] (30 



es - 3 P ) 



„ (1S -3P ), 

Z i/(A/=1) - 



rad MeV"2 



'(AJ=0) 

for nd. They yield two independent constraints on the five LECs. Absent any cancellation 
between different PV parameters we estimate that the rotation angles for both hydrogen 
and deuterium targets are of the same size, 



d0; 



PV 



dl 



[10~ 7 ■ ■ • 10" 6 ] 



rad 



m 



(VIII.l) 



There is no indication that the spin rotation observable is enhanced for a deuteron target, 
in agreement with Refs. jl^, 1^. Numerical stability analyses verify the theoretical findings 



of Ref. 



38( that no parity-violating 3NI is necessary at leading or next-to-leading order. 



The two calculations presented here join three others published in the same framework; 
^he longitudinal asymmetry in pp scattering [l| and two PV observables from the np -H- d'y 



system 
in Refs. 



2]. Note that other two-nucleon PV calculations using pionless EFT can be found 

1 n □ 

461 . 14 71 ] . To go beyond consistency to the realm of potential prediction, we 
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will next consider PV observables in Nd scattering, nd < — > 3 H_7 and pd < — > 3 He7 before 



considering heavier systems. 
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Appendix A: Partial- Wave Projectors 

In this Appendix we construct the projectors used to extract the desired partial wave 
state(s) from a given Nd s u state. The partial wave state is labelled [ 2S+1 Lj, I], where S 
is the total spin, L the orbital angular momentum, J the total angular momentum, and / 
the isospin. The projectors given here do not exhaust what is needed to do a partial wave 
decomposition for higher order or inelastic calculations; that would require constructing pro- 
jectors acting on NNN states as well. We employ the cluster-configuration basis introduced 



in Ref. 3jj, which also presents the S-wave projectors. 

The N aa field has two free indices: the SU(2) isospin index a and the SU(2) spin index 
a. The dibaryon field d\ is an isosinglet and spin-triplet ( 3 Si) with a free vector index i, 
while df is an isovector and spin-singlet ( 1 So) with a free isovector index A. To obtain a 
total 5 = 1 from an Nd s / t combination (indices suppressed), we can consider Nd s , which 
is already purely S = \ from the SU(2) decomposition | <g> = |. An S = ~ term is also 
available from Nd t via \ (g> 1 = § ffi \. The S = | piece is isolated by contracting with a 
Pauli matrix in spin space: 4^(o"j)^A raa (iJ. Summation of repeated indices is implied. The 
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normalization chosen will be explained at the end of this Section. 

An 5 = 1 term can only be obtained from Nd t . It has four degrees of freedom in spin 
space: one free vector index (the i in d\) with one free SU(2) spin index (the a in N aa ) gives 
six degrees of freedom, but the two that contribute to S = \ are removed by additional 
constraints on the projector. The most general form is given by 

Q i j = a5j + bo i Oj. (A.l) 

Requiring that the S = ~ and S = | projectors are orthonormal 

<Ti Q l 3 = = Q'j a j , Q\ Q k 3 = Q l 3 , (A.2) 

yields 

Q i 3 = S i j - = l\ 26 }- i e V ° k } » ( A - 3 ) 
where the second expression uses crVj = 5j + ie l j k a k . Orthogonality to the S = \ projector 

results in the required constraints that reduce the number of degrees of freedom to four. 
The projector onto S = | satisfies 

((/.,) Q J -- (A.l) 

For P-waves, the projector is obtained by combining L = 1 powers of the unit vector e 
in the direction of the d s /t center-of-mass momentum with the auxiliary field d s /t and the 
nucleon. In the case of the spin-1 field dt, this leads to total angular momentum 1 (8) 1 <gi | = 
§©§©!■ The projector onto J = I is: 



9 

Q 1j m ■= ^ 



(A.5) 



which is symmetric and traceless in (ij) and (kl) separately, 

SijQ ij k i = S kl Q l hi = , Q ij M = Q ji u , Q ij H = Q l3 lk , (A.6) 
{Q tj u) ] = Q kl tj , (A.7) 

and orthonormal to the S = | and S* = | projectors above: 

e>i Q^fei = Q lj ki a k = , Q m i Q lj ki = Q lj k i Q l m = , Q lj m n Q mn ki = Q lj u- (A. 8) 

Projectors onto definite states are obtained by contracting Oj and Q l j with d s / t N to extract 
the desired spin-state and multiplying by L powers of momentum. Any remaining free 
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indices are then contracted with the appropriate projector onto the desired total angular 
momentum, which is again given by aj, Q l j, or Q^kh an d with Kronecker-Deltas or e y as 
needed. For example, (ake h )(aid l t N) has the quantum numbers S = h,J 



I T 

2 ,1j 



1, i.e. 



2 Pi. Or, as a more complex example: 4 P3 mandates first coupling d\N to S = § using 
Q m i, resulting in Q m id l t N, and then multiplying by e^. The two free vector indices have 
to be contracted with another projector Q\, multiplied from the left, such that a J — | 
state results. This can only be achieved by multiplying with e nk m . The projector is thus 
proportional to Q l n e nk m ek Q m id\N . 

The same construction principle holds for isospin projections. States with I = \ can be 
generated by combining either d t (I = 0) or d s (I = 1) with the I = \ field N, since as 

| state, first contract all 



before <g> | 



\ and 1 <g> ~ 



| © \. To obtain an S 



spin- vector indices of d t N aa using <T/, or all isospin- vector indices of dfN aa using ta, leading 



to two equivalent terms 

5 a b 4 + (r A ) a b 5% df] N b ? = [a t d\ + r A df] N. 



(A.9) 



Here, all spin and isospin indices are listed explicitly in the first expression, while obvious 
index contractions are suppressed in the second. As two cluster-configurations exist, namely 



d t N and cLiV, it is convenient to follow Ref. 
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O = N, 



d} s \ 

u t,j> u s,B J 



( 



App. A.l] in decomposing all operators as 

\ bp 

(A.10) 



"7 

d-} 



N° 



0(Ndt -> Nd t )\ 0(Nd s -> Nd t ) J A 
\0(Nd t -> Nd s )f 0(Nd s -)• Nd s ) B A j 

Each operator is represented in the cluster-configuration basis by a 2x2-matrix which carries 
spin and isospin indices, and all operators act in the direct tensor product space spin ® 
isospin ® cluster. The operator in the cluster-configuration basis that projects onto the 



S 



component is 



{'Pd,iA) m ^aa 



1 

7! 



( n^ 

Oi 



mil 



(A.ll) 



On the other hand, only d t N contains an S — | 
define 

Q l j ( 




, / = \ component, so that it is useful to 



\l qj) act 



u 



mil 



Q^Vq withP ( 



q •" 



7 



(A.12) 
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Here, V q is the matrix projecting onto the only physical component of the S = | cluster- 
configuration matrix. The following relations hold: 

V* d A = (V dM ? = «V dM " , V dM V d A = l, 
. + (A.13) 
(V 3 V = V 1 ■ V k ■ V 3 ■ = V k ■ a V 1 ■ = 

\> q,ij ' q,j ' ' q,j 1 q,i ' q,i' d,lA ' qj u • 

The unit matrix in spin and isospin space is (l) b/3 a a — ^a^a- 

Below are the projectors as sources of the fields N and d s /t with the desired quantum 
numbers. First, in analogy with the auxiliary two-nucleon fields d s u, we introduce three- 
nucleon interpolating fields with the quantum numbers of the Nd s / t state. In symbolic 
notation, they are the cluster-configuration vectors representing sources T[ 2S+1 Lj, J] f with 
total spin S, orbital angular momentum L, total angular momentum J, isospin /, and 
appropriate vector and spinor indices {M}. The projection onto a definite partial wave in 
cluster-configuration space is then 

(T^L j7 Tf) m VM lA [ 2S+l Lj,I] (|) N . (A.14) 
The complete set of S- or P-wave projectors using auxiliary fields as source states is finally: 



P[ 2 Si,/ = ±] M = v d 



2 



I A 



V[ 2 Yi_J=\] lA = (a-e)V dMl 
V[ 2 Ps,I = \}\ A = V3Q\e k V dM 



P[ 4 P f ,J = i]^ = ^Q\ e nk m e k Q m l V q ^^[a i S k + a l 5 ik -Aa k 5 i l -5ie ik l ]e k V q , 
P[ 4 P f ,/ = i]^ = Q ij kl e k V q . 

(A.15) 

The spin-quartet projectors carry no isovector index; the normalization is discussed below. 

The two expressions for T^Ps] are equivalent since Q m \ contains two Pauli matrices and 
is multiplied with one more from the left, and products of Pauli matrices can be reduced 
to a sum of terms containing at most one Pauli matrix. The first form contains the spin- 
quartet projector P™ z = Q m {P q explicitly and hence is manifestly orthogonal to the spin- 
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doublet projectors. The second form contains the minimal number of linearly independent 
structures, derived by building the most general matrix A 1 m out of a 1 , <5] and e l ^ k . 

The orthonormalization condition was imposed as follows: Projectors to different partial 
waves are orthogonal after contraction over the auxiliary-field variable (IA) and integration 
over the solid angle element dfl e of the auxiliary field cm-momentum direction e: 

f dl e V[ 2S+1 Lj,I]^ lA (vf^L'j^I'A 1 * = V L^L'WS^S'VJ^J'VI^I'. 

An J \ J {N} 

(A.16) 

When the states are identical, integration and contraction must yield the identity element 
in the space of total angular momentum states. In the cluster-configuration basis notation: 

f ° r J= \ ihj ^eV[ 2S+1 L h J} lA ^[ 2S+1 Lr,I} lA = 1 , (A.17) 

t r O i for S — - 

for J = 6 -: -L dO e P[ 2S+1 L 3 _,I)\ A V^L^ = T' 2 ' (A.18) 

2 4W 2 2 [Q^V q for 5 = 1, 

for J = I: -j- f &l e V[ 2S+1 L,_,I] l i mn V^[ 2S+1 Ls , I] mn kl = Q i \ l V q . (A.19) 

1 47T J 2 2 

Appendix B: Assessing choices for the PC amplitudes of the class-II diagrams 

In this Appendix, we discuss in more detail the choices of the PC amplitudes t[X; k, q] of 
Eq. (1VI.13I) in the class-II contributions at NLO appearing in Eqs. ( IVI.10|MTT^VI.12p and 
shown in Fig.dJJ As discussed in the main text, the convolution kernel itself is already NLO, 
so that in principle a LO PC amplitude suffices. Figures [T5| [T6l and [T71 compare the depend- 
ence of the total amplitudes on the different treatments of t[X; k, q], benchmarked against 
the leading-order result and the choice of Eq. ( 1VI.13|) to include the 2 ' 4 S-wave amplitudes 
as NLO and the 2 ' 4 P-wave amplitudes as LO, see Eq. ( lVI.13p . 

Consider first the inclusion of the PC 2 Si-wave amplitude at LO instead of NLO (+ 
in Fig. [T5|) . while treating all other amplitudes according to Eq. ( 1VI.13I) . This leads to 
individual contributions diverging logarithmically and as g°- 23 - - ; as demonstrated on general 



grounds in Ref. [38| and numerically confirmed here. The problem is evident in the A 
dependence of the functions d and c of Eqs. flVI.19j) and (IVI.25j) . as shown in Figs. [T5l and 
[T71 This choice must thus be discarded as not properly renormalized. 

Next, it is consistent to include the PC 4 S3-wave amplitude at LO instead of NLO (x in 
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Figure 15: Cutoff dependence of the functions d[ 2 Si — > Y; coupling], Eq. ()VI.19j) . for dif- 
ferent choices of the order of the PC amplitudes in the class-II contributions of Fig. [TT] and 
Eqs. (1VI.10|ATTT]JVI12I) . ♦ (dashed lines): LO result, 3NI H from nd 2 Si-scatt. length; all 
other: NLO results with different treatments of PC amplitudes, with 3NI from 2 S±-scatt. length. 

2 

• (solid lines): PC amplitudes 2 Si as NLO, 2 ' 4 P as LO (Eq. (|V1.13|) 1: o (dotted lines): all as NLO; 
+ (dotted lines): all as LO. Linear extrapolations used to guide the eye. Notice the different scales 
on the vertical axes. Amplitudes which are identically zero are not displayed. Columns (left to 
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Figure 16: Cutoff dependence of the expansion coefficients d[ 4 S3 — > Y; coupling], Eq. (IVI.19I) . 

2 

for different choices of the order of the PC amplitudes in the class-II contributions of Fig. [IT] 
and Eqs. (|VL10^fVLTT^IVI.12p . ♦ (dashed lines): LO result; all other: NLO results with different 
treatments of PC amplitudes. • (solid lines): PC amplitudes S3 as NLO, 2,4 P as LO (Eq. (|VI.13|) ): 
o (dotted lines): all as NLO; x (dotted lines): all as LO. Linear extrapolations used to guide the 
eye. Notice the different scales on the vertical axes. Amplitudes which are identically zero are not 



displayed. Columns (left to right): partial waves S; 



2 P 3 , 4 Ss 



P3; rows (top to bottom): 



coefficients of Si, S2, T. 
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Figure 17: Cutoff dependence of the expansion coefficients c[g( x Y >] of the neutron spin-rotation 
angle in deuterium, Eq. (jVI.25|) . for different choices of the order of the PC amplitudes in the 
class-II contributions of Fig. Ql] and Eqs. (|VI10|rtVlTLjivT12p . Notation as in Figs. US and HH1 
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Fig. H~6]) . while treating all other amplitudes according to Eq. ( 1VI.13j) . This still determines 
the PV amplitudes at NLO since the kernel of the convolution already counts as NLO. With 
the exception of <i[ 4 S3 — > 2 P3 , S 2 ], however, the corresponding results are very close to those 

2 2 

obtained in a LO calculation of the PV amplitudes ( * ). Using the PC 4 S3-wave amplitude 
at NLO (• and o) induces a considerable shift of all functions <i[ 4 S3 — > Y; coupling]. As 
discussed in Sec. IVIB| PC LO and NLO 4 S| amplitudes predict values of the 4 S| scattering 
length that differ by about 25%, with the NLO result very close to the experimental value. 
The PV amplitudes, and therefore the <i[ 4 S| — > Y; coupling], are expected to be highly 
sensitive to the 4 S3 scattering length. This is supported by the plots in Fig.dHJ We therefore 
choose the NLO amplitude for the PC 4 S3 wave. Higher-order corrections from the 4 S3-wave 
amplitude are expected to be small, since the NLO expression for the scattering length is 
already in good agreement with experiment. 

Finally, it is consistent to include the PC 2 ' 4 P-wave amplitudes at NLO (o in Figures) 
instead of LO (• in Figures), while treating all other amplitudes according to Eq. ( 1VI.13f) . 
One might speculate that this leads to significant changes since the zero-energy effective- 
range parameter of the amplitudes, the scattering volume, is up to a factor of 2 bigger at 



NLO than at LO; see Ref. 



3l|. Instead, the amplitudes appear largely insensitive to the 



effective-range parameters of the P-waves. However, the spurious cutoff dependence already 
seen when the P-wave amplitudes are included at LO is widened and increased to a pole in 
the window A e [2000 . . . 4000] MeV in the total amplitudes. Its origin seems to be that 
the partially-resummed approach includes some contributions beyond NLO which need to 
be renormalized by 2,4 P-wave 3NIs which are not present at NLO. Outside that window, 
the functions c and d agree with the choice made in Eq. (1VI.13j) within the error-estimate 
of Sec. IVTEl 
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